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Abstract 

A class of continuous renormalization group flows with a dynamical adjust- 
ment of the propagator is introduced and studied theoretically for fermionic 
and bosonic quantum field theories. The adjustment allows to include self- 
energy effects nontrivially in the denominator of the propagator and to adapt 
the scale decomposition to a moving singularity, and hence to define flows 
of Fermi surfaces in a natural way. These flows require no counterterms, but 
the counterterms used in earlier treatments can be constructed using them. 
The influence of propagator adjustment on the strong-coupling behaviour 
of flows is examined for a simple example, and some conclusions about the 
strong coupling behaviour of renormalization group flows are drawn. 



1 Introduction 

The renormalization group (RG) method provides a way to attack, and in many 
cases solve, the infrared problem of models in field theory, statistical mechanics 
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and many-body theory. Infrared problems are ubiquitous and physically natural, 
since they are inherent to systems with gapless excitations, such as critical systems 
of statistical mechanics, massless field theories, and systems with spontaneously 
broken continuous symmetries. 

In systems of correlated fermions modelling metals, excitations around the 
Fermi surface are gapless. The attempt to treat interaction effects by naive per- 
turbation theory leads to a severe infrared problem in the form of divergences 
due to small denominators. The RG has been used for weak, short-range two- 
body interactions to cure this problem in perturbation theory and to classify the 
interaction terms by their power counting degree dlS. The one-body (or self- 
energy) terms, quadratic in the fermion fields, are the relevant ones, followed by 
the marginally relevant parts of the two-body interaction, which are quartic in 
the fermion fields. All higher terms are power counting irrelevant at weak cou- 
pling. The self-energy terms change the Fermi surface, which is the location of 
the singularity of the fermion propagator. Because the singularity has to be kept 
track of precisely, neglecting these terms leads to problems even in formal pertur- 
bation theory (provided it is done to higher than second order). For this reason, 
the self-energy terms are the most relevant ones. The marginally relevant parts of 
the two-body terms determine the quasiparticle interactions and give information 
about ordering tendencies and collective phenomena. The hierarchy of equations 
couples the evolution of the different terms. 

By the above arguments, a possible way to proceed is to assume that the Fermi 
surface given by the free propagator is already that of the interacting system. One 
can then take the dispersion function of a certain form, hence the interacting Fermi 
surface of a certain shape, assuming that this can be realized in an interacting 
system, and treat the remaining selfenergy effects by keeping them in the form 
of two-legged vertices. This also involves the assumption that there is no large 
correction to the quasiparticle weight. Both assumptions can, with considerable 
technical effort, be verified for the situation of smooth, positively curved Fermi 
surfaces (2l|3l[T0||. Technically, this is done in Q using counterterms to fix the 
Fermi surface. The counterterm function determines the shift from the Fermi sur- 
face of the free to that of the interacting system. The above realizability question 
takes the form of an inversion problem, which is solved under certain conditions 
in0. 

Although the importance of the self-energy terms is known and clearly ac- 
knowledged in the literature, a way to proceed in applied studies has been to focus 
on the study of the two-body interactions while neglecting self-energy terms in 
the flow. The main motivation for this is simplicity. Controlling the relevant self- 
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energy terms carefully is technically complicated already in the case of a smooth 
Fermi surface. The counterterm method of ||TJ provides a conceptually clear 
way of dealing with the problem, but it has been used only in a few applied stud- 
ies so far. 

In the two-dimensional Hubbard model, the case with Van Hove singularities 
is of interest for high-T c superconductivity [4]. For this case, a renormalized 
expansion using counterterms is developed, and some regularity properties of the 
self-energy are shown, in 0. However, the question of realizability of singular 
Fermi surface shapes has not been fully answered yet, although it is very important 
for deciding how generic Van Hove singularities really are in two-dimensional 
interacting systems Q. The effects of Fermi surface dynamics near Van Hove 
points can be rather nontrivial, since the Fermi surface is most easily moved in 
the vicinity of a zero of the gradient of the dispersion relation. This point was first 
clearly realized in 0, and it was shown that quantum fluctuations of the Fermi 
surface can arise and lead to non-Fermi liquid behaviour \7\. In this situation, it 
is useful to have a method in which the adjustment of the Fermi surface is done in 
the flow itself, because the concept of a non-fluctuating Fermi surface itself may 
make sense only down to a certain scale, below which the fluctuations emerge and 
may become dominant. 

The method discussed in this paper allows to use dynamically adjusting scale 
decompositions and to study Fermi surface flows. As a method, it applies more 
generally and can be used flexibly to adapt and adjust propagators in bosonic and 
fermionic (or mixed) systems, both for continuous and discrete flows. Variants 
of the method are used in mathematical studies of Fermi surface flows and Fermi 
liquid theory [3|. Also, the method has already been used to calculate the Fermi 
surface flow for the Hubbard model (see the Appendix of [ 8 ]). 

Various renormalization group schemes have already been used to study corre- 
lated fermion systems: Polchinski's original scheme [9|, the Wick ordered scheme 
ifTOirm . the 1PI scheme lfT2lfT3l l8l. Row equations for the 2PI functions were 
introduced iflU and studied for quasi- Id systems ffTSll . In this context one may 
well ask why one would want to investigate another scheme. In the following, I 
give some further motivations for this. 

As is well-known, all the above-mentioned schemes are equivalent on the 
level of untruncated hierarchies, because all of them contain the full information 
about all Green functions, as long as they are kept as infinite hierarchies. But they 
start to differ once the approximations that are necessary in nontrivial applications, 
are made, most importantly truncations of infinite hierarchies, gradient expansions 
of the generating functionals, or simplifying assumptions about the dependence 
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of the Green functions or vertices on momenta and frequencies. It is then useful 
to understand which features of a given scheme make it useful for a particular 
situation. 

All of the above-mentioned schemes include, in their untruncated form, all 
self-energy effects in the flow, albeit in a rather different way. In their original 
form, the Polchinski and Wick ordered schemes keep self-energy effects as ver- 
tices, while the 1PI scheme includes them in the most natural way, namely by 
having full propagators on the lines. On the other hand, the Wick ordered scheme 
can be transformed so as to have full propagators on the lines (see Section 4.7.2 of 
HH), and a similar transformation also applies to the Polchinski scheme. How- 
ever, the standard momentum space setup of the RG with a cutoff in momentum 
space near to the Fermi surface does not include an adaptive scale decomposition, 
and thus does not allow for the above-mentioned flow of the Fermi surface. This 
is the case in all schemes, so a correct procedure would require using countert- 
erms even in the 1PI scheme (in most applied studies, the Fermi surface shift was 
neglected, so that this problem did not appear). Alternatively, one can try not to 
use a cutoff on the fermion momentum ifTTl IT51 . but power counting has been 
done rigorously in more than one dimension only in the momentum space scheme 

CI [Ell- 
in two- or higher-dimensional applications, it is a serious practical problem 
to take the momentum and frequency dependence of the flowing vertex functions 
into account. The full propagator G of the 1PI scheme has support on all scales 
above the cutoff scale. This makes it necessary to extend the loop integrals over all 
of momentum space, which in turn requires having an accurate form of G, and of 
the higher vertex functions, over a wide range of momenta and frequencies. This 
has not been achieved so far. The situation in the other schemes is better in this 
respect: in the Polchinski equation, all internal lines carry a derivative of the cutoff 
function. Therefore, at scale e, the propagators get evaluated only at values of the 
momentum k with \e(k)\ ~ e. The integration procedure for Polchinski's equation 
makes the equation nonlocal in the scale parameter in the standard way lfl"9ll. but 
a significant advantage is that whenever a momentum region |e(A;)| ~ e' further 
away from the Fermi surface contributes to a loop integral in the flow equation, it 
does so only via a propagator on scale near to e' , which was already calculated in 
an earlier stage of the flow, and which does not need to be adjusted in every further 
integration step. In the Wick ordered scheme, all propagators are supported at or 
below the cutoff, so only small neighbourhoods of the Fermi surface occur in the 
loop integrals at low scales. Thus the propagators and vertices are needed only 
close to the Fermi surface and at small frequencies. 
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In studies and talks by the Gottingen group lEUll . it was pointed out that in 
simple models, the 1PI scheme has a better strong coupling behaviour than the 
other schemes. Namely, the 1PI flow equation is accurate in these toy models at 
much larger coupling values than one might believe at first, far beyond the weak 
coupling regime, and it remains accurate in a region of intermediate to strong 
couplings where all other schemes fail. The dynamically adjusting scheme is ex- 
amined for this toy model in Section 4. It turns out that propagator adjustment, 
namely putting self energy effects into the denominator instead of keeping them as 
vertices, is a necessary, but not sufficient, ingredient in stabilizing the equations 
at strong coupling. The Polchinski scheme still fails at strong coupling, while the 
dynamically adjusted Wick ordered scheme is accurate up to arbitrarily large cou- 
plings. The analysis also clarifies the reasons for the success of the 1PI scheme. 
In Section 4, 1 also give a more general discussion about strong initial couplings, 
and show that for reasons that are simple, but specific to the RG method, the 
"weak-coupling" RG flow is a well-defined starting point at strong coupling also 
in "real" models (such as the two-dimensional Hubbard model), contrary to what 
one may believe at first, and to what is stated in many places. Whether this can be 
used efficiently in computations depends on certain relative rates of growth; this 
is under investigation and will be discussed at the end. 

2 Theory 
2.1 Setup 

The partition function of a general bosonic or fermionic theory with source fields 
^ (the generating function for the Green functions) is 

Here the fields Q x are indexed by elements X of a finite set X. X is taken finite to 
make the generating function mathematically well-defined. In many applications, 
X = T x I, where T is a finite space-time lattice and / is a finite set of internal 
indices, but this assumption will not be needed for the general theory. If X comes 
from such a lattice regularization of a continuum theory, the limit of lattice spacing 
to zero and volume to infinity can be taken later in the RG equations (and after a 
volume factor has been divided out in the field-independent term). Readers who 
trust that all this can be done properly (or who don't care) can also think of a 
continuous space T right away. 
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For bosons, the &x are real-valued; for fermions, each <d x is one of the gener- 
ators of a Grassmann algebra. Spinors and vector fields are included in this setting 
by an appropriate choice of /, and complex bosonic fields as well, by regarding 
them as two-component real fields (with the corresponding index again part of /). 
The commutation relations are 

&X&Y = (&Y&X (2) 

with ( = 1 for bosons and ( = — 1 for fermions. The integration measure is VQ = 
Ylxex d0x, an d O) = a Ylix&x ^ x®x with a > 0; (when considering a 
lattice approximation, a is the volume of the unit cell of the lattice). The action 
of the operator A is (A Q) X = a ^2 y (A ) X y®y', the matrix elements (A ) x ,y 
are in general complex. An operator C is called (^-symmetric if 

(*,ce) = c(c*,e), (3) 

(i.e. Cx,y = (Cy,x)- hi the bosonic case, Re C is also required to have nonnega- 
tive eigenvalues, so that for all 6, 

(6, ( Re C)0) > (4) 

(in the case of complex fields viewed as two-component real fields, this corre- 
sponds to the usual condition that the hermitian part of the covariance is positive 
definite). Assume that the A in (JTJ is (-symmetric. If A is invertible, both 
A and the propagator C = A^ 1 are (-symmetric. The normalization factor M 
is chosen such that Po(0) = 1 for V = 0. For bosons, I also assume that V 
is bounded below and that ~(0, AqQ) + V(0) grows faster than linearly in all 
1 0x| for large \Qx\- Under these hypotheses, the integral in (JTJ is a convergent 
finite-dimensional integral in the case of bosons. In the case of fermions, it is a 
linear functional on a finite-dimensional space. The theory can also contain both 
bosonic and fermionic fields. Every formal functional integral can be regularized 
so that it takes this form. 

The logarithm of P generates the connected Green functions of the theory. 
Assume that A is invertible and let C = Aq 1 . A shift in the integration variable 
gives (using the symmetry of the quadratic part under 6 — > —0) 

P (tf) = e ^*>*) J d^ c {Q)e- v ^ +c ^ (5) 

with dyUc*(0) the normalized Gaussian measure with covariance C. In this form, 
one can take a limit where some of the eigenvalues of C go to zero; in the bosonic 
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case this corresponds to restricting the integral to the subspace ker C. Typically, 
this situation arises by introducing cutoff functions that vanish strictly in parts of 
momentum space; see the following examples. 

Although the equations derived in this paper work for the general case set up 
above, it is useful to keep two prototypical examples in mind. 

The first example is scalar field theory on a finite lattice T C Z 3 where |/| = 
1, (A ) is the discrete Laplacian with (say) Dirichlet boundary conditions, and 
Vo(0) = J2xer V{®x), with V a polynomial that is bounded below. In that case, 
(A )x,y = a>o(X — Y) with Fourier transform a (k) = 2^] 3 =1 (1 — cos/cj) for 
k 7^ 0. The momentum k = is removed by the Dirichlet condition; of course, 
other boundary conditions can be used, too, provided the Laplacian is restricted 
to the complement of its kernel. For details, see, e.g., Chapter 2 of ifTol . 

The second example is the regularized partition function for spin 1/2 fermions 
on a finite sublattice A of Z d , d > 1. Here Y = {0, |, 2|, . . . , (n - x A, 
I — {1,— 1} x {1, — 1}, with the first index being a charge index and the second 
the spin index. For X — (c, £) G I x V, with £ = (a, r, x) and X' = (c', £'), with 

(A) - ( Ao &^ \ (ti\ 

{Ao)x > x '-{-M?,o o) cc , (6) 

where A (£, £') = S a:Ct /a (T — t',x — x') and the Fourier transform of a is of the 
form a (u, k) = iuj — e (k), where the Matsubara frequency uj takes values in odd 
multiples of tc times the temperature T = and eo is a real-valued function that 
has a nontrivial zero level set, the Fermi surface So — {k : eo(k) = 0}. Vo( l I / ) is 
an even element of the Grassmann algebra with Vo (0) = 0. In the simplest case, Vo 
is a quartic polynomial in the fields, e.g. Xi v q(X 1 , . . . , X A ) ^ x ^ ■ ■ ■ ^x 4 > 

but it can be much more general. The number n is an auxiliary quantity coming 
from an application of the Lie-Trotter product formula. To recover a representa- 
tion of the partition function of the many-fermion model, as given by a trace over 
Fock space, one must take n — > oo first. For details, see Chapter 4 of IfTol . 



2.2 Renormalization group 

The renormalization group comes into play when the inverse of A becomes un- 
bounded, or even fails to exist, in the limit \X\ — > oo, or when the inverse remains 
bounded, but its norm becomes nonuniform in an essential physical parameter. In 
the above examples, the Fourier transform diagonalizes A and one can read off 
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that the eigenvalues go to zero in the limit of interest. In Example 1 above, taking 
T to be a cube of sidelength L, k { = riilix / L with n« £ N ; clearly, a (k) goes to 
zero for infinitely many n = (ni, n 2 ,n 3 ) in that limit, so that C = clq 1 becomes 
singular. Similarly, in Example 2, one gets zero eigenvalues when the temperature 
T — > 0, and when fc £ So- At small positive T,C — a 1 is of order 1 /T. In both 
cases, C, and thus C itself, fails to be square integrable in the limit T — > 0, so that 
integrals over a square of a propagator (and higher powers of it) are ill-defined. 
This is not just mathematical pedantry but physically relevant - if in the case of 
the many-fermion system, C were square integrable at T = 0, there would be no 
superconductivity, because the superconducting gap is driven by the log (3 from 
J \C\ 2 . Similarly, the singularity of 1 in the bosonic case leads to anomalous 
decay exponents of correlation functions. Hence in the following, a guideline for 
getting a useful RG equation in these systems is that squares of unregularized 
propagators are to be avoided. 

A standard approach to obtain a renormalization group equation ffOl in the 
above examples is as follows. Introduce a scale parameter s > 0, an s-dependent 
energy, say e s = e e _s , which goes to zero in the limit s — > oo. Choose a fixed 
decreasing C°° cutoff function x< with x<( x ) — 1 for s < 1/4 and x<( x ) — 
for x > 1, and set x> = 1 — X<> so mat X< + X> — 1 i s a smooth partition of 
unity. Decompose the propagator into an infrared part D s and a part supported at 
higher scales, C s , via 

C(k) = C s (k) + D s (k) = X > (^^) C(k) + X < (^^) C(k). (7) 
and define the effective action 

= -log J d/i a (e)e- v ( 0+$ ). (8) 
The partition function is then expressed in terms of the effective action as 

P (*) = e ^ c *M j d m (6)e-^( e+ *). (9) 

The effective action Q s is a function of a fields which, when expanded in the fields, 
has a quadratic term. In infrared problems, this term is relevant in the RG sense. 
That is, the order of the singularity may get changed, producing the anomalous 
exponents, in the scalar field theory example, or, in the many-fermion system, 
the changes to the quadratic term also shift the location of the singularity of the 
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propagator, the Fermi surface, in Fourier space. The above scale decomposition 
does not take this into account, because it zooms in on the fixed Fermi surface 
5*o. Therefore, with such a fixed choice of scale decomposition, counterterms are 
needed to avoid divergences [2]. It is possible to justify the use of counterterms 
by an inversion theorem [2], but it seems desirable to have a simple method for 
moving the quadratic terms from Q s into the propagator, in a way that the shift 
of the singularity gets included in the flow automatically. The flow introduced in 
this paper allows to put self-energy terms into the propagator and hence to de- 
fine a convenient adaptive scale decomposition both in continous (and in discrete) 
renormalization group flows. It is such that an analogue of continues to hold, 
but with a propagator that changes dynamically in the flow, and Q(s) replaced by 
an interaction from which all, or part of, the quadratic part is removed. 

Finding a scheme that implements an adjusting scale decomposition in a dif- 
ferential equation is not completely straightforward. On graphical grounds, one 
expects the full propagator to be G = (A — E) -1 , where E is the Dyson self- 
energy. Thus the first idea may be to make the ansatz 

c W = x> (4f£M) m 

with A(s) = A — E(s), E(s) to be determined. C(s) is well-defined because x> 
is nonzero only on a subspace where A — E(s) can be inverted, and the cutoff 
function now adapts dynamically because it depends on E(s). The propagator 
of the unintegrated fields must be defined as = C — C(s) to ensure that 

P(s, \&) = Po(^). But this has a singularity at the noninteracting Fermi surface 
for any s, as well as one developing on the new one, so (apart from being a very 
unnatural object) it does not vanish for s — > oo, so that not all degrees of freedom 
get integrated out. Then this flow would not give useful information about the 
model, even if it converged for s — > oo. Trying instead to keep the form of @, 
but putting 

D(s) = AW - x< p^M) 

leads to disaster: The derivative contains a term — A(s)~ l 9 ^^ A(s)~ 1 x<, 
which has a square of the inverse, but no infrared cutoff because x<( x ) — 1 f° r 
x < 1/4. In all cases where the self-energy shifts the Fermi surface, ^ on 
the singularity set, so that effectively, a square of the propagator appears. Because 
this propagator is not square integrable, this scheme leads to divergent integrals 
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already in the RG differential equation at fixed s. The method developed in the 
following sections avoids these problems. 

2.3 Conditions on the adaptive flow 

I now return to the general field theoretic situation of Section |2~T1 and consider the 
generating function in the form Q, i.e. 

P (tf) = e 3(A.*,w) J df i Do (Q)e- v ^ e+D ^ (12) 

Instead of D = C = A 1 1 make the slightly more general choice 

D = A \ (13) 

where xo can be chosen as an operator that cuts off very large energies, chosen 
such that D is (-symmetric, and for bosons, nonnegative (but xo — 1 is also 
allowed). In applications, such an ultraviolet cutoff allows to start the flow at a 
scale where degrees of freedom with very high energies have already been inte- 
grated over. The flow is set up by posing restrictions on the general, s-dependent 
partition function 

P(s, *) = e K(s)H(G(s)^) J ^ (s)(e)e -v( s ,e+s (sW _ (14) 

Here K, G, V(-, S, and D to be determined as functions of s. For fixed s, 
K(s) is a complex number, S is an operator, G and D are (-symmetric operators, 
and V(s, \l/) is a function of the fields which, for bosons, has to be such that 
the integral exists (see below). 

I impose the following general conditions. 

HI K(0) = 0, G(0) = D(0) = S(0) = D , and V(0, ip) = V (ip). 
H2 For all s > and all P( s , *) = P(0, 
H3 D(s) -»• as s -»• oo. 
These conditions have the following simple interpretation. 

1. P(0, vp) = Po( v P) is the generating function of the model of interest. 
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2. The generating function P, and hence all correlation functions, are invariant 
under the RG. 

3. In the end, all degrees of freedom are integrated out. 

The field-independent term —K(s) is (up to a factor equal to the free en- 
ergy. The term involving G(s) naturally arises once one starts adding terms to 
the denominator of the propagator. This is most easily seen in the case when one 
requires that V(s, should contain no quadratic term in Then, if the limit 
s — > oo exists and D(s) ^ in that limit (as required by condition 3), 

P(oo, = e*°° + 5(G<»*.*) e -v^(Soo*)_ (15) 

Because V(s) has no quadratic term for any s, the same holds for Vqo = lim V(s). 

Therefore the term involving Goo is the only quadratic term and hence Goo is 
the full propagator. The limit exists e.g. if the temperature is above all critical 
temperatures [1701 |2"T1 |3) or if the Fermi surface is asymmetric under reflections 

ED. 

Thus, the full propagator G is built up gradually in the RG flow. For finite 
s, (TRjl has the usual RG interpretation that D(s) is the propagator of particles 
with interaction V(s), but now V(s) can be chosen not to contain any self-energy 
terms. Thus, even if the limit s — > oo cannot be taken because the flow is accurate 
only up to some so (see, e.g. JUJ0D), the effective propagator and interaction can 
still be used as an input for the theory of the excitations with energy below e SQ . 

The above conditions do not determine K(s) and the functions G(s), V(s, 
S(s), and D(s) uniquely. They are just the minimal conditions one wants to 
impose. The freedom of choosing the various functions can be used to impose 
the constraint that the quadratic part of V(s) should vanish, or that a certain part 
of it be removed from V(s) and put into D(s) instead. 

The flow (J9j) satisfies H1-H3, but it restricts to the propagator D(s) = XsC 
with the fixed scale function Xs of (J7J), which does not adjust to a shifting location 
of the singularity. This condition is relaxed when one drops © in favour of ()14[). 
In particular, one can choose the functions such that the propagator D(s) incor- 
porates all, or part of, the selfenergy corrections. In applications, this means that 
one can directly study the flow of the Fermi surface, the flow of superconducting 
gaps and the like, without having to use counterterms. 
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2.4 The RG differential equations 

I seek only solutions to the conditions H1-H3 for which K, S, G, D and V are 
differentiable in s. Then HI and H2 are equivalent to the initial-value problem 

§- s P(s, V ) = (16) 
K(0) = 0, G(0) = £>(0) = S(0) = A), V(0,V) = W), 

from which the RG differential equation will be obtained. 

For a (-symmetric operator A, the Laplacian in field space is defined as 

A A = (5v,A 8v) (17) 

where 5^ = a -1 ^. 

Theorem 2.1 Let s \— > Q(s) &e given, with Q(s) (symmetric for all s and 
bounded for all s. Assume that there is a solution (A, x) to the initial-value 
problem 

A{s)=- X (s)Q{s), A(0)=A , x(0)=xo, (18) 
such that s h- > x( s ) is differentiable and the low-energy propagator 

D{s) = A{s)- l x(s) (19) 
and the fluctuation propagator 

F(s)=A(s)- 1 -^x(s) (20) 

exist and are (-symmetric, and that D(s) is nonnegative for bosons. Assume that 
K and V satisfy the initial-value problem 

j>( a ,¥) = e v ^ l -A F(s) e~ v ^ + l -(Q( s )^ ^) + K - C -Tr {DQ), 
V(0,tt) = Vo(tt), #(0) = (21) 

(Y/zms in particular V(s, is differentiable in s and at least twice differentiable 
in ^) and the condition that V(s, 0) = Ofor all s. Set 

S(s) = A(s)~ 1 xo (22) 

and let G be the solution to 

G(s) = S(s) T Q(s)S(s), G(0) = D . (23) 

Then the functions K, G, D, S and V solve the initial value problem {73}, hence 
they satisfy HI and H2. 
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This theorem is proven in the next subsection. 
Remarks. 

1 . D must be (-symmetric, and for bosons it must be nonnegative. A way to a- 
chieve this is to restrict to (-symmetric A(s) that commutes with its adjoint 
A*(s), and to take as a function of A*A(s), so that the two commute, 

e.g. x( s ) = X< ( ^^T^^ \ ■ The conditions on A are satisfied in the above 
examples and in most other interesting applications. 

2. The dependence of the cutoff function x(s) on A(s) allows for an s-depen- 
dent, adaptive scale decomposition. This makes showing the existence of 
a solution to (fTSJ) nontrivial. This existence problem will be discussed in 
more detail in Section l2~9l 

3. The fluctuation propagator F(s) given in (|20|) is now indeed proportional 
to the derivative of the scale function, thus has an infrared cutoff. Therefore 
the problem mentioned after (|TT|) is absent from this flow. 

4. The equation for the full propagator G can be rewritten as 



Thus in general, G(s) ^ A(s) 1 , also in the limit as s — > oo. This must be 
so for general reasons, as will be explained in Section l2~6l 

5. The strategy to satisfy H3 as well is to choose a solution for which the 
support of x(s) becomes empty as s — > oo, and for which the singularity 
of A(s)^ 1 as s — > oo is such that D(s) = A(s)~ 1 x(s) still vanishes in the 
limit. 

6. Theorem 12 . 1 1 gives a family of solutions to the conditions HI and H2 para- 
metrized by Q. Moreover, it will become clear in the proof that even given 
a fixed Q, there are many other possibilities to choose the equation for V, 
and as explained, there is a further freedom of choosing \ as a function 
of A when seeking a solution A) to (JT8J). Thus, there seems to be an 
enormous freedom in solving H1-H2. However, not every choice has good 
analytical properties, and most solutions diverge already at a small value of 
s. To have a bounded solution on a maximal interval for s requires a careful 




i 




-i 



(24) 
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choice of Q(s). For instance, in the Fermi surface problem in d > 2, Q 
must take the Fermi surface shift into account fully; whether it also needs to 
contain corrections to the Fermi velocity and the coefficient of cu depends on 
finer details, such as the shape of the Fermi surface. If this is done properly, 
the limit s — > oo of the flow exists at temperatures above the mean-field 
transition temperature [3|, and H3 is satisfied. 

7. Compared to the Polchinski equation obtained from (J5J) by differentiating 
with respect to s, which reads 

g s = e G ^A D e- g % (25) 

(with D s = Cx< as in (J7j)), there are two changes. First, there are terms 
of order and 2 in the fields ^ on the right hand side of (|2T|). which can 
be chosen to cancel the field-independent and quadratic part arising from 
the first term on the right hand side, so that V(s) starts with terms of order 
4 in ^. Second, the fluctuation propagator F(s) now includes selfenergy 
corrections from the integration up to scale s (if Q(s) is chosen to absorb 
only selected parts of the quadratic part generated by the integration, e.g. 
only cancelling the Fermi surface shift, then the propagator A(s), and hence 
F(s), contain only those). 

Thus fl2Tj) has a simple and natural interpretation: the propagator is allowed 
to change dynamically, and the selfenergy (or part of it) can be put into the 
denominator. Although the RHS of (J25J) has the same structure as the first 
term on the RHS of (j2T|) . the change of the quadratic part, and the ensuing 
replacement of C s by F(s) leads to a major change in the behaviour of the 
RG equations. Some aspects will be discussed in Section|4| 

2.5 Proof of Theorem 2.1 

To make the proof more readable, I first give an algebraic derivation which applies 
both to bosons and fermions, and in a subsequent section provide the mathematical 
justifications of all steps. 

2.5.1 Derivation of the RG equations 

Define the convolution with the Gaussian measure as 

(/i D * = / dfi D (e)F(e + $) (26) 
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The partition function is 

P(s, rj) = e K(s)+ms)v,v) ( ms) * e-vW) (5(5)77) 
If is differentiable in s and vf, 



-W(s,S(s) V ) 



dW 
ds 



with 



#(s) = 5(s) S(s)- 1 and V H = (HV, 6*) . 



The equation 4^ = implies 



with 

Using 

gives 



Q(s) = {Sisy'f G(s)S(s)-\ 



II — -^^HD+DHT 



(27) 



(28) 



(29) 



\K{s) + ^(Q(s)% *) + P^w + J^J * e- v «) (*) = (30) 



(31) 
(32) 



(k(s) + \(Q(s)V, *) + + ^ e^ A ^e- v ( s '*' = (33) 
Acting with e v ( s >*)e~3 A £>« n (g2J and using 



(34) 



e -3 A »i(Q¥,¥)e3 A " = ~(Q*, *) - | Tr (DQ) + ^A DQD - P OQ (35) 
(recall ( = 1 for bosons and ( = — 1 for fermions) I get 

V = e v (k - I Tr (£>Q) + i(Q¥, *) + ^A F + 2^ e~ y (36) 



with 



F = D-HD- DH T + DQD, B — H — DQ. 



(37) 
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To simplify the equations, I remove the dilation term V B , i.e. choose S so that 

H = S^S = DQ (38) 
This implies F = D - DQD, so that 

VW = K-^Tr(DQ) + ^^) + e v ^A^ DQD e- v ^ (39) 

This equation makes it natural to regard Q(s) as an input, and to seek a solution 
for D,V, and K. Once this is done, the remaining functions are obtained by 
integrating the equations S = DQ S and G = S T QS. 

Let D be of the form dHJ) with A(s) and x( s ) me solution to dHJ). Dif- 
ferentiation with respect to s implies that F(s) is indeed given by (}20|) . With 
the choice (fT~9|) . the amputation operator S can now also be determined: S = 
DQS = -A A^ 1 S, with initial condition S(0) = A(0)~ 1 Xo has the unique solu- 
tion S(s) = A(sy 1 xo stated in (|22|) . 

Thus Theorem l2.1l is proven algebraically for bosons and fermions. Adapting 
the derivation to the case of a theory involving bosons and fermions is straight- 
forward, provided that the interaction V is an even element of the Grassmann 
algebra. 

2.5.2 Mathematical remarks 

In the above derivation, it was assumed 

1 . that V and P are indeed differentiable in s and in \P 

2. that (J32J) holds, and 

3. that the commutator equations hold, when applied to e~ v . 

Item 1 follows from the existence of the solution to the initial value problem for K 
and V, as stated in the hypotheses of Theorem l2.1l (and already explained there). 
Concerning item 2, note that for any polynomial II in the fields, 

(fintf * n) (*) = e^ A «wn($) (40) 

was proven in ifTol both for bosons and fermions, but e~ v is in general not a 
polynomial. Similarly, the commutator equations fj3*4*)l and (J33J) hold when applied 
to polynomials in ^I^, but again, the above proof involves an application to e~ v . 
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Consider first the case of fermions. In the present setup, the Grassmann al- 
gebra is finite-dimensional, so every function of the Grassmann variables is a 
polynomial. Thus V(s, 6) and e ±v ^' e ^ are polynomials in 6 and P(s, ^) is a 
polynomial, hence infinitely differentiable in Thus (J32|) holds, and (}34|) and 
P5|) hold as well. Moreover \ is differentiable in s by hypothesis of Theorem 12. 11 
and A is differentiable in s because it solves the differential equation (fT8|) . so D 
is differentiable in s as well. By P2"j) . the convolution in fb^) is a polynomial in 
D(s) and in V(s, •). Thus by expanding V in homogeneous parts of degree 2m in 
\& (see also Section l2~6l) . one easily proves by induction on m that V, and hence 
also P, is differentiable in s. 

For bosons, e~ v ^ s '^ is not a polynomial, but, under our assumptions, it has a 
fast decay for |\&| — > oo. This decay and some analyticity could be used to justify 
first truncating e _v to a polynomial and then taking a limit. A different, easier 
way to adapt the argument is to avoid (J32j), i.e. go back to (}30|) and use integration 
by parts-arguments of the type 



* / d(M D (* ~ ©) ^(©) = / <W* - 6) - F(0) (41) 



where J-(s, \&) is the difference between the left and right hand sides of f|39|) . 
Keeping D(s) > in the derivation and obtaining the general case of nonnegative 
D by a limit, an equation of this form implies that 6) = 0. This yields the 
same differential equation as the algebraic method of the previous section and thus 
completes the derivation for bosons. 

In conclusion, all steps in the derivation of the RG equations in the previous 
subsection are justified mathematically. In the case of fermions, the smoothness 
of V(s, in s and \& follows from easy general arguments. For bosons, it is part 
of showing that the solution to the equations for V and K exists. 

2.6 Structure of the dynRG hierarchy 

For brevity I call the RG differential equation for V, K and A, \ derived above 
the dynamical RG or dynRG equation. Expansion of V in powers of the fields 
gives a hierarchy of equations for the m-point functions. In this section, I assume 




to bring the equation into the form 




(42) 
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that V is even and choose Q and K such that V does not contain any quadratic 
and constant part, which gives the hierarchy with the simplest graphical structure. 
The expansion in the fields 

V(tf) = V (m) (^) (43) 

m>2 

with iX" 1 * 1 homogeneous of degree 2m in ty, leads to the hierarchy (shown graph- 
ically for Q, and V® in Figure[TJ 



m-i 1 i 
yM = £ _(^y (n) , F <5*V (m+1 - n) ) - -A F V (m+1) , m > 2 (44) 

n=2 

(g(s)*,*) = A FW V (2) (s,*) (45) 
K = C - Tr (DQ) (46) 

By dHJ), (03) reads k = -§ Tr = § Tr log A, so it integrates to 

= -~ ( Tr log A(s) - Tr log A(0)) (47) 

For fermions, £ = — 1 and Tr log A = Pf A, where Pf (A) denotes the Pfaffian of 
the antisymmetric matrix A. For covariances of the form @, Pf (A) = (det A) 2 , 
so that 

e*« = detA(s) (48) 
det A(0) 

For bosons, ( = 1 and Tr log A = det A. Thus for real bosons, 

-1/2 

.det A(0) 7 

For complex bosons with a hermitian covariance there is a factor 2 which squares 
this to det A(0) / det A(s). 

Thus the relation of the operator A to the free energy density is simple: if 
A(s) — > A for s — > oo, the free energy is (modulo a conventional factor of the 
inverse temperature) simply given by 

K = -| ( Tr log A - Tr log A ) (50) 
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e ^ )=( detA( g )V (49) 



This also clarifies why the equation (J23j) for G cannot simply be G = A~ l : the 
exact relation between the full propagator G and the free energy involves the 
Luttinger-Ward functional. It is a nice feature of the dynRG flow that the relation 
between A and the free energy K is so simple. 

Clearly, the diagrammatic structure of the hierarchy (see Figure[l]) differs from 
the Polchinski hierarchy obtained from (|25|) by the absence of one-particle func- 
tions (or two-point functions) which get absorbed into A via the equation (ff5|) for 
QandJIU). 

The standard iterative solution of this hierarchy produces a graphical expan- 
sion. It should be noted that from a purely graphical perspective, the dynamical 
propagator adjustment does not make any of the functions one-irreducible. This 
is obvious for the six-point and higher functions, but also holds for the graphs 
contributing to the four-point function V*- 2 * 1 and to Q: taking the tree diagram 
of the six-point function and joining two legs of one of the vertices to a tadpole 
produces a one-reducible contribution to taking two legs of the other vertex 
and joining them to a tadpole gives a one-reducible contribution to Q. The point 
here is that these two-point insertions (the single-scale insertions defined in 0) 
are harmless because they do not generate singularities. In the Fermi surface ex- 
ample, the size of such an insertion at scale s is by power counting (see [ 10 1) of 
order e s ~ e~ s , while the size of the propagators on the adjoining lines is of order 
e^ 1 ~ e s . Thus the combination of two propagators and such a two-point insertion 
is of order e~ s e 2s = e s , which is the same size as a single propagator. However, 
if one inserted a two-point graph that is integrated from to s, its value would be 
~ Jq ds'e~ s = 0(1). Then the combination of this insertion with two propaga- 
tors is of size e 2s , hence e s larger than it should be. It is contributions of this type 
that lead to a problem with unrenormalized expansions. Here, this problem does 
not arise because renormalization is implemented by a changing propagator. In 
summary, for the analytic questions, namely which terms are large and which are 
small, the reducible graphs contributing to V' 2 -* and Q are inessential. 

2.7 Wick ordering and self-consistency equations 

Expanding in Wick ordered polynomials provides essential simplifications in the 
mathematical analysis of the RG equations ifTOl IT6l : moreover the Wick ordered 
scheme has been useful for calculations in two-dimensional systems ffTTI l23l . A 
distinctive feature of the Wick ordered equation is that all propagators are sup- 
ported on scales below e s , so that only low-energy degrees of freedom enter the 
integration. In the many-fermion example, this means that all momentum inte- 
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Figure 1: The first three equations of the dynRG hierarchy. The first equation 
determines Q, which enters the right hand side of the differential equation for the 
inverse propagator. The dots denote derivatives with respect to the scale parameter 
s. All internal lines carry a fluctuation propagator F. 

grals are restricted to a neighbourhood of the Fermi surface. Therefore projec- 
tions to the Fermi surface, which are needed to solve the equations numerically, 
are better under control than in other schemes. An further advantage is that struc- 
tures of overlapping loops El [KB |24l| become explicit in the equations and make 
the study of improvements to power counting much easier. These improvements 
are needed for showing Fermi liquid behaviour in higher dimensions and control- 
ling two-particle interactions. 

A natural (but not the only) way to fix the Wick ordering covariance is to 
require that tadpoles be removed. In the case of the fixed scale decomposition (J7J), 
this is easy to do : the condition to cancel tadpoles is to have -j^R(s) = —C s , 
because then the linear term in the RGDE drops out [ 1 1 . In this case, the solution 
to that equation is simply R(s) = D s , with the D s of (J7J). In the adaptive scheme 
pT|) . the corresponding equation is 

^R(s) = -F(s) (51) 

which is not solved by R(s) = D(s). The reason for this is simply that a Wick or- 
dering that removes all tadpoles requires information about the propagator at later 
scales s' > s, which requires solving self-consistency equations. In theoretical 
studies, this self-consistency can be solved by fixed-point arguments. Alterna- 
tively, one can choose the Wick ordering covariance such that it leads to the above 
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improvements without cancelling tadpoles exactly. 
The Wick ordered equation is obtained from (J2TJ 



*(.,»)= - iA FW V( S ,*) + I(^),^ S ^*) 



+ ±(Q(s)%*) + k-^Tr(DQ), (52) 



by setting 
Then 



V(s, = e-5 A «M W(s, *). (53) 

1 

2 



V(s, = - A A(s) V(s, *) + e -* A *W VV(s, *) (54) 



and the resulting equation is 

W(*,*0 = -Ia f{s) _ r{s) W(s^) 

+ [e' A ^)^A^W(s,*i)W(a,* 2 )] (55) 
+ ~(Q(s)*, v]/) + £( a ) - | Tr [(!>(*) - fl( S ))Q(a)]. 

Here the usual combinatorial device E51 of introducing two copies ^ and of 
the field has been used to rewrite the quadratic term, and [ . . . ] denotes evaluation 
at \t>! = \tr 2 = ty. 

The choice of R as the solution to 

R(s) = F(s) (56) 

makes the linear term in the equation drop out, as in ifTOl . There are some further 
conditions on R that were automatically satisfied in the case without propagator 
adjustment in ifTUIl . The support of R needs to go to zero as s — > oo, and R needs 
to satisfy certain power counting bounds to guarantee that Wick ordering does not 
introduce singularities. For this reason, rewriting the solution in terms of an initial 
condition as 

R(s) = R(0)+ f 'dt A{t)- l x{t) (57) 
Jo 

is not useful. Instead, requiring R(s) — > for s — > oo, the solution is conveniently 
written as 

dt Aity'xit) (58) 
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However, in contrast to the situation without adjustment of propagators (where 
R(s) = D s ), this is really a self-consistency equation because it requires knowl- 
edge about the propagator at later RG times. For general purposes, such as order- 
by-order arguments and power counting, it suffices to know general properties of 
R, and to solve for R inductively while solving the RG equations. In practical 
applications, one needs to use approximate solutions R instead of R, for which 
R — R is less singular than R itself. 

2.7.1 The hierarchy of Wick ordered functions 

The expansion 

W(s,tf) = W (m) (s,*) ( 59 ) 

m>2 

with W ( ' m - ) homogeneous of order 2m in ^ gives the following hierarchy. For 

m > 2, 

yvM( Sj tf) = -^A F _ R W {m+1) (s,y) 

+ E E ^ [(Ag ,2 V Ag^> W^*, ^ 
for m—1 

+ EEx K A n 2) ) tA F 2) W M (*>*i) W^>( S ,* 2 )] (61) 

£>2 M=2 

and for m = 

AT = | Tr [(£> - R)Q] 

+ EE^T K^R A ) e ^F A W< 1+ ^( S ,* 2 )] (62) 

^>3 m=2 

This hierarchy is shown graphically in Figure |2j The number of loops of the 
corresponding graphs is i. As is explicit in the above equation, the contributions 
quadratic in W start at two-loop for Q and at three-loop for K. The choice 
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Figure 2: The first three equations of the Wick ordered dynRG hierarchy, with the 
choice R = F and the truncation that = for m > 4. In the truncation 

where also = 0, only the first of the three terms remains on the right hand 
side in each equation. 

F = R removes the term that is linear in W. A graphical analysis of the iteration 
implies that with this choice, there are no one-reducible contributions to and 
to Q. Moreover, if the one-loop term is dropped from the eq. for the four-point 
function, then all graphs contributing to and Q are 2PI. 

2.8 Ladder summations and symmetry breaking 

If the Wick ordering covariance R is chosen as the solution to R = F that vanishes 
for s — > oo, as above, and if the hierarchy is truncated by setting VyX" 1 - 1 = for all 
m > 3, the resulting equation for becomes 

= \ I {U^Stf) WW(«,* a )l (63) 

The structure of this equation is analogous to that obtained by a similar truncation 
from the 1PI hierarchy after the 1PI hierarchy has been rearranged in Katanin's 
KWi way: Eq. (25)) in |E71 is to quadratic order in the interaction of the same form 
as (J63j) . It is also analogous to equations (4.176) - (4.178) in [ 1 6 1 because there a 
product of propagators D S D S appears in the flow equation as well; however, there, 
D s does not contain self-energy corrections). Note that in contrast to the 1PI case, 
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the structure of the scale derivative of a square of the Laplacian (which is essential 
for the solution of the selfconsistency) appears in the Wick ordered dynRG in a 
natural way; no rearrangement of the hierarchy is necessary. 

It is thus clear that, when restricted to a particular channel, the equation is 
easily solved by a geometric series in analogy to the procedure in Section 4.5.4. 
of iTToTl . In the present framework, restricting to a channel means that (possibly 
after a change of variables), the equation for the coefficient function Wx,y,x',y' of 
has the structure 

W X ,Y,X',Y>(s) = ^ Wx,Y,U,v{s)U-U,V,U',V'{s)W U ' > v' >X ',Y'{ s ) ( 64 ) 

uy,u',v 

which has the structure W = WJ1W with the product defined in the way evident 
from ({54*)) . i.e. regarding pairs (U,V) as matrix indices. Because ({"""")) contains 
a derivative of the square of the Laplacian, II is indeed explicitly given as the 
derivative of a matrix IT. The solution of ([64)) is 

w(s) = w(o) (i - [n(s) - n(o)] w(ofj . (65) 

That integrating the RG equation leads to such a resummation was, together with 
the self-energy equation, crucial for accessing the symmetry-broken regime with 
a fermionic flow in *""71 . The above structure suggests that this can also be done 
using the dynamic Wick ordered hierarchy developed here. 

In this example, one resums in the particle-particle channel, so II(s) is simply 

ftu,v,u',v'(s) = R(s)u,v R(s)u',v (66) 

where R(s) denotes the non-charge-invariant propagator with A(s) given by the 
matrix in (37), (38) of l"2"71 . multiplied by the infrared cutoff function. However, 
there is a crucial difference to the 1PI scheme, which I now discuss at the example 
of the flow of superconductor done in [27|. By definition, the Wick ordering 
propagators R(s) are supported at small scales below e e~ s , so Hindoo R(s) = 0. 
Thus lim^oo n(s) = 0, and ({55)) gives for s — > oo 

= W(0) [1 + U(0)W(0)]- 1 . (67) 

That is, the solution of the channel RG equation (j65)) does not lead to the final, 
but the initial value of the propagator in the solution of the four-point function. 
Hence, in the solution of the BCS model given by this truncation, it is the initial 
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gap that appears in the solution for the four-point function. This seeming contra- 
diction (when viewed from the perspective of ||2"71 . where the initial gap is put in 
as a symmetry-breaking parameter that is subsequently sent to zero) is resolved 
when one remembers that in the Wick ordering scheme one also has to Wick order 
the initial interaction. The requirement that the initial interaction is Wick ordered 
but has no quadratic part in the fields leads to the Hartree-Fock self-consistency 
equation. In the reduced BCS model, this equation yields the exact solution for 
the superconducting gap. That is, in the Wick ordered dynRG, the initial Wick 
ordering of the bare interaction already gives the exact solution of the model, so 
the gap is constant and equal to the initial one. The flow of then gives the 
exact vertex. The two-loop contributions to Q, as well as the higher contributions 
to , all vanish in the thermodynamic limit, as explained in the Appendix of 



2.9 Relation to discrete RG flows 

In this section, I discuss a discrete analogue of the flow differential equation, 
and show that the differential flow with a dynamically changing propagator can 
be obtained as a limit of this discrete equation. The successive changes of the 
propagator in the discrete scheme take the form of a continued fraction expansion 
for operators. 

This approach gives a specific way to check the existence of the solution to 
the equation (fTSj). which is fundamental for the flow developed here. A full math- 
ematical proof of existence is rather nontrivial, but the following gives a clear 
prescription how to proceed and what to check. Discrete flows are interesting by 
themselves since they are used in the mathematical studies in this field. 

2.9.1 Adaptive iteration 

Consider again the effective action Q for a system of fermions with propagator D 
and interaction V, defined by 



G2. 




(68) 



Q is of the form 



Q(F, V, *) = K(F, V) + -(Q(F, V)V, *) + I^(F, V)(V), (69) 
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where /C is independent of ^, Q defines a quadratic form that gives the quadratic 
term of G in \1>, and Z^- 4 ) contains all the higher powers of the fields. 

By successively splitting the covariance in two pieces, integrating over fluc- 
tuations, and shifting quadratic parts into the measure, I will now recursively 
construct a sequence of propagators and interactions (D n , V n ) n >i with the prop- 
erty that the quadratic part of V n vanishes. Assume that Vo has no constant and 
quadratic parts. Split D — E + F , then 

e 6(A,,Vo,*) = J d ^(^ x ) e G(F ,Vo,*+*)_ (70) 

By Q(F , V ) splits into a constant term K = 1C(F , Vo), a quadratic term 
Qo = Q(F ,V ), and a term Vi = J(- 4 )(£(F , V )). The Gaussian measure 
is dfi Eo (^) = AfiE^e'^^^D^, with the normalization factor M for the 
Gaussian measure being a power of determinant or Pfaffian of E , so the quadratic 
term given by Q can be shifted into the measure, and get 

d/W*) e^«°*) = Af(l - E Q ) d^*) (71) 
with the propagator 

D 1 - 1 = E - 1 -Q , (72) 

which now includes the corrections to the self energy given by Q . Note that Q 
is not the proper selfenergy S but instead the connected amputated two-point 
function. The relation to the selfenergy S coming from the integration over \l/ 
is 

g = (l-S Fo)- 1 So (73) 

If one takes the limit of a strict support condition for F , Q = S outside the 
support of F , i.e. at low energy scales. 

In writing f)72|l. I have assumed that the RHS is invertible. This condition 
will be further discussed below. An obvious variant of this formula holds if one 
decides to put only part of Q (e.g. the part that shifts the Fermi surface) into the 
propagator. 

In the convolution integral on the right hand side of (|70)L appears instead 
of so that 

e -G?(A),Vb,*) _ gfif+^(*,Oo*) f dfi Dl (^) e (*-Qo*)+Vi(*+*) ( 74 ) 
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withfTo = K+\ogJ\f(l — E Q ). Completing the square in the Gaussian measure 
and using the identity 1 + DxQ = DxEq^ 1 , gives 

G(E + F Q ,V ,$) =Ko + l^,(Qo + QoD 1 Q )^)+g(D 1 ,V 1 ,D 1 E - 1 ^). (75) 
Iteration of this identity by splitting D x — E x + Fi and proceeding as above gives 

n— 1 n—1 

g(D ,v ,$) = J2^ + 2^ Sl ^' Ql ^ 1 + Dl+lQl)Sl ^ 

1=0 1=0 

+ g(D n ,V n ,S n $) (76) 

with recursively determined propagators D n , interactions V n and amputation op- 
erators 

n-1 

S n = Y[D k+1 E-\ (77) 

k=0 

Consequently, 

log P(D , V , v)=K n + ^(G nV , rj) + G(D n , V n , S n D oV ) (78) 

with 

n-1 

G n = D - YsiSiDofQ^l + A+iQO^A). (79) 
i=o 

If the splitting in E n and F n is chosen such that D n — > as n — > oo, and if lf n , 
V n , S'n, and G n converge to limits V, 5, and G in that limit, 

logP(A), V ,v) = K + i(G?77, 77) + V(^). (80) 

If at each iteration step the full Q n was put into the measure to define D n+1 , G is 
the full propagator of the model, and V generates the connected m-point functions 
of the model with m > 6, with amputation given by S. The interpretation at a 
finite n is similar, except that there still remains an effective interaction for the 
degrees of freedom that have not yet been integrated over. 

In the following, I discuss two specific prescriptions for the splitting, which 
cover most of the applications, and derive the corresponding formulas for S n , D n , 
and G„ . 



27 



2.9.2 Difference cutoff 



Assume that the original covariance D is given in terms of an invertible A as 
D = Aq 1 , and that A commutes with its adjoint A* . The operator \A \ 2 = A* A 
is hermitian and positive and therefore has a spectrum that is a subset of IRq = 
{x E H. : x > 0}. For n > 1, let Xn + Xn = 1 be a C°°-partition of unity 
on Wq, where Xn is decreasing and Xn is increasing and both are strictly positive 
functions. We also require that 

^T^<1, (81) 

Xn{x) 

the idea being that the Xn provide lower and lower cutoffs as n increases. A 
possible choice is Xn(x) = (1 + e 7 ™^ -1 )) -1 and % n = 1 — Xn, but the details of 
this choice do not matter for the moment. The operator x n (A*A), defined via the 
spectral representation of | A\ 2 = A* A, is a positive operator that cuts off the parts 
of the spectrum outside the support of Xn- It is not a projection because we chose 
a continuous partition of unity instead of a step function. With this choice of x as 
pointwise strictly positive, Xn(|^4| 2 ) is even invertible. 

Given Aq and the interaction V , the sequence is constructed as follows. Set 

E = A Q - l xii\A Q \ 2 ). (82) 

Then F = D -E = A - 1 Xl (\A \ 2 ). Calculate Q = Q(F ,V ) and V 1 = 
J^ 4 )(F ,V ), and define 

A 1 =A -xi(\A \ 2 )Q . (83) 

If A 1 is invertible, set A = A^ 1 xi{\A \ 2 ) . Then £>f 1 = E^ 1 - Q . In general, 
for n > 1, and given V n and an invertible A n , set 

D n = A n - 1 x n (\A n - 1 \ 2 ), E n = A n - 1 Xn+i(\A n \ 2 ), (84) 

and F n = D n - E n . Determine Q n = Q(F n , V n ) and V n+1 = l^(F n , V n ), and 
set 

A n+1 = A n - Xn+ii\A n \ 2 ) Q n . (85) 

If A n+ i is invertible, the iteration can be continued; otherwise, the iteration stops 
at n. The thus defined D n and E n satisfy 

AT 1 = E~\ - Q n _ u D n+1 E~ l = A~X x A n , (86) 
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so that 

N-l 

H D n+l E- 1 = A N l A Q . (87) 

n=0 

The condition that A n is invertible has to be checked in every step because it is, 
in general, nontrivial, and its failure can mean that there is some instability. We 
discuss this in more detail in the application to the Fermi surface deformation 
below. 

There is another, less obvious, property of this scheme: in the fluctuation prop- 
agator 

Fn = (xn(K_x| 2 ) - Xn+l(K| 2 )), (88) 

the difference of cutoff functions does not necessarily give a positive operator. In 
a bosonic theory, this would already make the fluctuation integral ill-defined. For 
fermions, lack of positivity is not a problem for the definition of the fluctuation 
integral, but numerically, positivity is important. Roughly speaking, the difference 
of cutoff functions becomes negative if the change of the propagator is so strong 
that one starts integrating backwards. 

One can modify the cutoff function by replacing A by PA, where P is some 
projection, whenever it appears as the argument of a cutoff function. For instance, 
in the fermion models, one can keep the cutoff function independent of the Mat- 
subara frequency variables. 

2.9.3 Positive cutoff function 

The following alternative scheme is more straightforward but does not admit a 
continuum limit if one wants to keep a continuous partition of unity. 

Given a single partition of unity x + X — h with x decreasing, and a strictly 
decreasing sequence of positive numbers (e n )„>i, let Xi = X (''4f _ ) an( * Xi — 
Set 

E = A^xi, F = A 1 Xl . (89) 

Then E + F = D , and both have positive cutoff operators. Again, calculate 
Qo = Q(F ,V ) and V 1 = I^(F Q , V ), and set A 1 = A - X iQo- If A 1 is 
invertible, set D\ = A^ 1 . 

We now proceed in exactly the same way, by splitting D\ = D± (x + x) — 
Ei + F\ and so on. This builds up a product of cutoff functions, which remains 
positive because each factor is positive. More precisely, for n > 1, and given 
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Aq,..., An, let 
and set 

D n = A~ l x„, E n = A- 1 Xn + i- (91) 

Then both E n and 

F n = D n -E n = A-'xn X ( (92) 

V e n+l / 

have a positive cutoff factor. After calculation of Q n , set 

A n+X = A n - Xn+lQn (93) 

and, if A n+1 is invertible (which is the case if Q n is small enough), D n+X = 
A'^Xn+i- Again, this sequence satisfies (|55jl. 

2.9.4 The summed iteration 

Both of the above schemes can be described as follows. We have a sequence of 
invertible operators A n and positive operators x n , such that A = A, xo = h and 
D n = A~ l Xn, E n = A~ x x n+ i, and A n+l = A n - x n+ iQ„, with Q n given by the 
fluctuation integral with propagator F n = D n — E n and interaction V n . Iteration 
gives 

n— 1 n— 1 

£(A>,V ,$) = + (Q, + Q,A +1 Q,)^r%^) 

+ £(D n ,V n ,,4;%$). (94) 

This equation has an easy interpretation: the Q on the RHS contains a fluctuation 
integral over fields with propagator D n and interaction V„, and the external field 
is amputated with A n instead of Aq, as it should be, since D n is A~ l times a cutoff 
function. The constant and quadratic terms are the contribution to the free energy 
density and the full propagator. 

The generating functional is thus given by 

log P(A), V , rj) = K n + i(G n 77, V) + G(D n , V n , A^v) (95) 
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where K n is the field-independent term and 

n-l 

G n = D - AT\Qi + Qi D i+iQi) A 7 l (96) 

1=0 

is the contribution to the full propagator up to step n in the RG iteration. 



2.9.5 Limit of a differential equation 



I now use the "difference cutoff' scheme of Section 12 .9 .21 to obtain a continuum 
limit of the RG sequence. Let eo > be an initial energy scale, let e n = e^—ne and 
choose Xn(x) = x(-t) where x + x = 1 is a fixed partition of unity on M.t with 
strictly positive functions x and x- Since there is a one-to-one relation between 
n and e n and we want to take a continuum limit, we now label the sequence by 
e„ instead of n, and denote A n = A 6n , V n = V 6n , etc. Since in the limit e — > 
nothing is integrated over, Q en is of order e: 



Qe n =eQ' tn . 

Because e n ±i = e n =F e, the fluctuation propagator at e 



is, by 



A: 



X 



\A 



- X 



\A f 



e J \\ e ~ £ 
Thus, up to terms that vanish as e — >• 0, F' e = jF e is given by 

d_ 

We 



(97) 



(98) 



* = A" 1 -x' ( 



(99) 



Note that the e-derivative acts only on the cutoff function, but not on the inverse 
in front. The differential equation for Q e and the other functions is obtained by 
doing the fluctuation integral to 0(e). Using (|52"|). recalling that F n = F' e, and 
dropping terms of order e 2 , 1 get 



Q(F n ,V n ,<S>) = log ((1 + A F Je v *W) 



l og (e v "(*)(l + e- v "^A Fen e v "W)) 



Thus 



0(ei^,V e ,$) =V e ($)+e 



Af'H($) + 



5$ ' 6 5$ 



(100) 



(101) 
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hence the increment is simply given by the right hand side of Polchinski's equa- 
tion. The main difference here is, however, that V e has no quadratic part. Writing 
V e = Ve^ + V e , extracting the quadratic part of Q, and taking the limit e — > 0, 
gives the differential flow equation as follows. The change of the inverse propa- 
gator is given by 

JU e = Q'M 1 -^ 1 ) (102) 

where Q' e is determined by 

i($,g' e $)=A F ,V( 4 )($) (103) 

The interaction part obeys 

= A F ,Vp 6) ($) + \ (<J*V e , F'MV e ) (104) 

These are the equations obtained directly in the continuum in Section 12.41 (with 
the change of variables e — > s = log e /e). 

2.9.6 Wick ordering in the discrete scheme 

The discrete scheme also provides insights about the choice of Wick ordering in 
the limit of a differential equation, and it is of independent interest. The point 
where things start to differ from the non-Wick-ordered scheme is the definition 
of D\. Instead of (J72j) . D\ is now defined such that the integral in (|74|) is rewritten 
in terms of a Wick ordered interaction. More precisely, first set $ = in (f74|) and 
require that 

/ ^- Qo) -i(*> vm = ^ J d^M^ imm (io5) 

where = e~^ AD F denotes Wick ordering of F, as in fTTOl . and 

V { - 2) Wi = 0, (106) 

i.e. Wi(\l/) does not have constant or quadratic parts in \E'. This last condition 
leads to a self-consistency equation for D 1 whose solution differs from (j72]) . By 
definition of the Gaussian measure, this means that 

logA/^- 1 - Qo)- 1 ) - \{ (E, 1 - Q )V, ¥) + Vi(¥) 

= logjV(Di) +Ki- -{D^, *) + e^Sw^) (107) 



32 



Wi is obtained by using 1 = V^- 2 ^ + T^- 4 ) to rewrite 

V 1 = e"^ (V^ 2) + 7^ 4) ) e^i V x . (108) 
Since V^- 4 ^ projects out all terms constant or quadratic in the fields, this gives 

Wi = P { - 4) e5 A ^i Vi 

K x = p(°) e ^ Ac i Vi - ^A Dl V {2) e^ A ^ Vi (109) 

and 

D" 1 = Eq" 1 - Qo + P (2) e5 A °i Vi (1 10) 

where for a function of the fields F ((/>), A = P^Fis the operator in the quadratic 
form given by V (2) F, i.e. V {2) F{§) = \(A $, $). Equation (fTTOl) is the self- 
consistency equation for D\. It takes the form of a generalized Hartree-Fock-type 
equation because the action of ea A °i is to create self-contractions of the vertex 
given by Vi with a D\ propagator, and this feeds into D^ 1 in the way that is 
standard in self-consistency equations of mean-field type. Now that Wi and Di 
are given, one obtains the rewriting of the integral in ff74]) for $ ^ by standard 
shifting of Gaussians and redefinitions of the amputation term that acts on the 
external field $. 

3 Fermi Surface Flows 

In this section, I briefly discuss the Fermi surface flow in d > 2 in the scheme 
without Wick ordering. A discrete Fermi surface flow with Wick ordering is de- 
fined and controlled mathematically for d = 2 in [3j. This is one of the examples 
that motivated the method. Although different methods, such as the T-flow of 
ifPTl have been used to calculate RG flows, and other regularizations, such as a 
cutoff on the frequencies only, are possible, the momentum space flow remains the 
only one where precise mathematical estimates for power counting of all Green 
functions, not just the first few truncations, have been given Q- 

Let the covariance and interaction be those of the standard many-fermion 
system with short-range interaction V , as described in the Introduction, with 
n — * oo. (for more details, see |[16|, Chapter 4). In Fourier space, the covariance 
is a multiplication operator, namely the inverse of 

A)(<?o,<?) = i<7o - e Q (q), (111) 
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times the unit matrix in spin space. Here e is the dispersion relation for free 
particles, given by the hopping amplitudes of the model. By standard (formal) 
diagrammatics, the full propagator is the inverse of ig — £o(q) — cr (9o> q), with a 
the fermion self-energy. One natural choice for would then be to take A as the 
full propagator, with a understood to depend on s (but without a cutoff function). 
In the following, I shall discuss another choice, where only the frequency zero 
part of the selfenergy is put into the denominator. This then implies that V still 
has a quadratic part, but one which vanishes at zero frequency. In cases where the 
field renormalization is a marginally relevant parameter, one can also put in the 
linear part of a in q . Thus the freedom in chosing Q and A allows to adapt to the 
situation at hand, while keeping the denominator as simple as possible. 
So, let 

A(q ,q) =\q -e(s,q). (112) 

In terms of the scale-dependent self-energy cr s (q , q), e(s, q) = e (q) + a s (0, q), 
so that the difference 

°s(qo,q) = o- s (q ,q) - (r s (o,q) (H3) 

vanishes at q = 0. Let x be a decreasing cutoff function with x(x) = 1 for 

x < 1/2 and x(x) — for x > 1. Set 

(another choice would be to use a cutoff on the frequencies as well: x( s > ko, k) = 

^ , fcg+e(s,fc) 2 



With this choice, V(s) still has a quadratic part, which is of the form 

fflc$(k)«i(a,k)^(k) (115) 



with 

v 1 (s,(0,k))=0 (116) 

(here k = (ko, k), and J dk — (3' 1 J2 kg( z Mf f with Mp the set of fermionic 
Matsubara frequencies). Thus with this choice of moving only frequency-inde- 
pendent parts of the selfenergy into the denominator, there remain some quadratic 
terms in the action. But by construction, these terms vanish on the Fermi surface 
(see (11161) ). and hence do not cause divergences. 
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The choice (jl 14|) leads to an adaptive scale decomposition because x is sup- 
ported in the neighbourhood 

{k: \e(s,k)\ <e s } (117) 

of the set {k : e(s, k) = e (k) + cr s (0, k) = 0} which one may regard as a flow- 
ing Fermi surface. In the standard, nonadaptive scale decomposition, the cutoff 
function is taken as Xo(s, k) = x(eo(p) 2 e~ 2 ), whose support shrinks to the free 
Fermi surface as s — ► oo. Consequently, the nonadaptive scale decomposition can 
be used only if a counterterm is included EJ HI - or if the Fermi surface shift is 
ignored, as is sometimes the case in approximate studies. 
In the present case, the terms in (J43j) are 

/m 

n=l 

V m (s\k 1 ,...,k m ,k' lJ ... J k / m _ 1 ) (118) 
Equation ([T8|) can now be rewritten as 

e(s,p) = Q(s,p)x(s,p) (119) 

with 

Q0,p) = J cfkF(k) V 2 (s | (0,p),k;k) (120) 
The fluctuation propagator 

F(k) = A(k)- 1 ^ X (,,fc) 



A(k)" 1 X' ( 4(e(a, fc) 2 + e(a, fc)e(«, fc)) (121) 



e 2 7 e 2 

s / s 



contains a term proportional to e, because of the adaptive nature of the scale de- 
composition given by this choice of x- 

To discuss the mathematical properties of the resulting equations, it is useful 
to introduce the "tadpole" operator. For a function k) it is the linear operator 

(Uf)(p) = 4 fdq Hp; q)x' (^4^-) -■ M {ll2) 

e 2 J \ e 2 / i?o - e(s, q) 

acting on functions of q. With 

$ v (p,k)=V( S |(0,p),k;k) (123) 
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and 



$ v ,e(p,k)=V(s|(0,p),k;k) e(k) 



(124) 



fjl 19|) becomes 



(1 - £<& v , e ) (e) = >C$ v (e 2 ). 



(125) 



The operator £ depends nonlinearly on e, but it is linear in e. 

Thus already the differential equation for the change of e involves an operator 
inversion , due to the additional term in the fluctuation propagator in (jl21j) . hence 
coming from the adapting scale decomposition. The solution strategby involves 
the following steps. 

1. Show that 1 — £$ Ve has a bounded inverse on some open set of V and e. 
This is essentially a restriction on the size of V 2 — it will fail if V 2 gets so 
large that its backreaction on e becomes of order 1. A sufficent condition is 
\\C\\ < 1 in a suitable norm. This allows to rewrite (|125|) as 



Because C depends on e, it is necessary to find a set of e on which ||£|| < 
r\ < 1 uniformly. The existence of this set already requires conditions on 
eo; the set is then essentially a ball around eg. 

2. The general theory of ordinary differential equations then implies existence 
and uniqueness of the solution for all s in some bounded interval, provided 
that F(e) is Lipschitz. Proving that this condition is satisfied is again non- 
trivial, in fact it is closely related to questions of possible symmetry break- 
ing. 

3. To get global existence (i.e. for all s, one needs to prove that the flow actu- 
ally never leaves the ball around e . 

4. In steps 1-3, V was still regarded as fixed; in the true system of equations, 
V 2 depends on s and the RHS of the equation for V 2 depends on e as well. 
Thus one needs to show that V 2 stays bounded. 

There is a similar set of steps for the discrete Fermi surface flow. Steps 1-3 are 
already nontrivial, but the hard part is including 4, because it requires dealing with 
the hierarchy as a whole. For the differential equation, an all-order proof that steps 
1-4 work can be given along the lines of ifTUll . A nonperturbative proof, i.e. with- 
out truncation on the hierarchy, remains an open problem, mainly because to this 



e = F(e) = (l-U v J- L U v (e 2 ). 



(126) 
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day there are no good nonperturbative bounds for the continous RG equation ||2"81| . 
which would be needed in step 4. Up to now, this step requires an integration over 
short scale intervals and corresponding estimates on discrete RG transformations 
E9l . Such a proof is given in the Wick ordered approach, for two-dimensional 
Fermi systems, in [3 1. In particular, the boundedness of V2 necessary in step 4 can 
be shown in the temperature range of the Fermi liquid criterion of ifTUl . 

These mathematical points also correspond to physically natural problems. 
The question of the existence of the inverse of 1 — C is, in physical terms, basic to 
the consistency of the method itself: the RG flow is parametrized by energy scales 
on the (changing) kinetic part of the action, i.e. the part quadratic in the fields. 
The nonquadratic interaction part is not taken into account when labelling the 
scales. However, it influences the change e of the quadratic part, and the equations 
become inconsistent when this change becomes larger than the maximal kinetic 
energy scale itself. Note that it is this condition that is important in the flow; the 
typical interaction energy, measured by the coupling constant, can be much larger 
than the scale e s , and this situation indeed occurs at low scales in all nontrivial 
models because e s — ► for s — ► 00. It is a nontrivial fact of fermionic field theory 
that the method can be controlled by convergent expansions even in that situation 
J29). The invertibility of A(s) at all scales s is related to this; the s-dependent 
Fermi surface must never leave the region where x<(e,7 2 e(,s, A;) 2 ) 7^ 0, otherwise 
an inconsistency manifests itself as divergences in the equation. This condition 
can be satisfied because power counting implies that the selfenergy at scale s is 
of order e s times the coupling constant (in fact, improved power counting [TU) 
gives a better s-behaviour of the bound, which is needed to control derivatvies of 
the selfenergy in the flow). 

The counterterms used in [2| can be obtained from the adaptive flow: if e(k) = 
lim^oo e(s, k) exists, the difference e — e can be taken as the counterterm K of 
0. This can be used to avoid solving the inversion equation 0. However, the 
question whether there is a one-to-one relation between the free and the inter- 
acting dispersion function and Fermi surface, which is answered by the inversion 
theorem, is not solved this way. 

The need to invert 1 — C for a tadpole operator C also arises in the 1PI scheme 
when one reexpresses the single-scale propagator S in terms of the full propagator 
G in the equation for the self-energy E 112611271 . With the tadpole operator defined 
in a natural way, the equation for the selfenergy in the 1PI hierarchy becomes 

t = £ r{2) (S) = £ rC2) (G - G±G) (127) 

where I^ 2 ) is the 1 -irreducible two-particle vertex and G is the full Green func- 
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tion. This equation can be rewritten as 



(1 - £ r( 4) GG )(£) - C T (4)(G). 



(128) 



More generally, the occurrence of such operators can be expected in any scheme 
where the propagator changes, be it with an adaptive scale decomposition, as done 
here, or without (as in the usual 1PI hierarchy). 

The flow discussed here can be adapted to a flow where the density, not the 
chemical potential /i, is kept fixed. This variant was used in the Appendix of [ 8 1 
to calculate the flow of the Fermi surface at fixed density in the Hubbard model. 

4 Strong coupling behaviour in a toy model 

In this section, I study the flow equations for a very simple example, namely the 
integral over one complex variable, corresponding to the partition function of a 
zero-dimensional complex bosonic field, 



with a > 0, Re A > (for simplicity, I take A real in the calculations below.), 
and h a complex source field. For tp = u + if, dip A dip = 2idudt> . Thus it is 
simply a two-dimensional integral. 

This analysis is motivated by ll20ll . where it was observed that the 1PI RG 
equations give an accurate solution for a function similar to Po(h) even for large 
A . This is unexpected at first sight because the diagrammatic approach to these 
equations may suggest that they are useful only for small A . Moreover, it is in- 
teresting because only the combination of RG and irreducibility seemed to induce 
the good behaviour at large couplings, namely (1) in contrast to the 1PI scheme, 
the other schemes studied in ll20ll . namely Polchinski's original hierarchy and the 
Wick ordered hierarchy of ifTUll . were accurate only at small A , as one would ex- 
pect from small-coupling schemes and (2) perturbation theory for the 1PI vertices 
T m fails as badly as in the other schemes when the coupling constant is not small. 
Thus it is not merely irreducibility, but a different feature of the 1PI RG hierarchy, 
that makes for the difference. 

A closer look at the equations reveals a first crucial feature: the 1PI equation 
contains the selfenergy o in the denominator. When A gets large, a becomes of 

1/2 

order A after a very short flow time already, and (because the signs conspire 




(129) 
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well in this toy model) the factor A in the numerator is balanced by a factor A 
in the denominator, which prevents the right hand side of the flow equation from 
getting very large. This suggests that a dynamic adjustment of the propagator, 
such as discussed in this paper, improves the behaviour of the other hierarchies as 
well. This is the case, but the above explanation is not the full story, as we shall 
see in the following. 



4.1 Asymptotics of P at large Ao 

Although we have written P as a function of h, it is understood that it also de- 
pends on h. Indeed, by the U(l) symmetry of the measure it depends only on 
hh = \h\ 2 . A Gaussian transformation gives the representation 

P (h) =a [ -£= e~ r2 / 4Xo (130) 

from which analyticity in \h\ is obvious. Write 

P o (/0 = ^7r m (a o ,Ao)W 2m (131) 

m>0 

Rescaling the integration variable in ()129|) . one gets n m (q , An) = A 2 F m ( ^2= 
where F m (£) is analytic in £ for |£| < 1. Explicitly, 

This convergent expansion gives the behaviour at large A as 



P (0) = vroK Ao) = ^ + 0(f) (133) 

* vAn An 



and 

p o( h ) e «u|A| 2 +«j2[fc|*+o([fc| B ) ^34) 



^o(O) 

with 

wi = 1/2 + 0( ^ W2 = " (h ~ s 1 A ° 1 + 0(A ° 3/2) - (135) 



Thus for large A , the 2m th moments of the integral, which play the role of the 
2m-point functions in this toy model, become small. In the following I discuss to 
what extent this can be reproduced by the RG with propagator adjustment. 
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4.2 The dynamical RG hierarchy for the toy model 

The usual shift gives 

P (h) = e^/ a ° * e-) (136) 

with the notations 

d ^ ) = ^ a ° e " a ° M2 ' (/i*F)(0) = |dMy?)F(0-^ (137) 
and w{tp) = A |y?| 4 . As in the general case, write 

P(s, h) = e fc W + ff«l fc l a (/i, (s) * e- w W) (138) 

and impose k(0) = 0, g(0) = a^ 1 , d = Oq 1 , er(0) = Oq 1 , v(0,(p) = w((p), as 
well as dP/ds = 0. The RG equations are derived using straightforward differ- 
entiation, in analogy to the general case. Expand 

v(s, V ) = Y,v m (s)\ V \ 2m (139) 

m>0 

then the equations read 

v = k-dq-fv u vi = q + f(vl - 4v 2 ) (140) 

and for m > 2 

v m = f ^ m + 1 ~~ A*) VWi-*i -(m+ l) 2 v m+1 j (141) 

The initial interaction u> has no constant and quadratic part. We now impose that 
the function v has the same properties, by requiring v = and v\ = 0. This gives 
equations for q and k: 

q = Afv 2 and k = dq. (142) 
Finally, we fix the flow by setting 

d(s) = a( S )- 1 X (s) (143) 

where a(s) > is to be determined and x( s ) > a given function of s that 
decreases from x(0) = 1 to = (here we take s = 1 as the final value of the 
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flow parameter. In this simple example, we do not need to make x depend on a). 
We require a and x to be nonnegative to have the same for d, so that the Gaussian 
measure is really a measure. Now set 

a = - X q (144) 

Then / = a _1 x- Note that / is negative because x is decreasing; the sign of / is, 
however, unimportant because it never appears in a measure, only in a Laplacian. 
With our initial condition, v 2 starts out positive; v 2 decreases and a increases in 
the flow. 

The system of equations reads 

k = a- 2 xx^v 2 (145) 

a = ~XX^~ 1 v 2 (146) 

v 2 = -9a^xv3 (147) 

v 3 = a^x (Avl - 16t» 4 ) (148) 

v 4 = a- 1 X (12v 2 v 3 -25v 5 ) (149) 

tm-l \ 
^ v(m + 1 - fi)v IJ ,v m+1 -p - (m + l) 2 v m+ i J (150) 



The summation in ()150|) now excludes \l — 1 and // = m because Vi — 0. The 
equation for a can be rewritten as 

§- s (a(s) 2 ) = -Av 2 (s) §- s (x(s) 2 ). (151) 

The "amputation" factor a(s) is now determined by the analogue of condition 
(J38j): putting aa^ 1 — dq = and inserting ()143|) and ()144|) gives o v b = —a~ 1 d, 
thus, by the initial condition, a(s) = a(s)^ 1 . Moreover, a glance at ()145|) and 
f|146|) reveals that k = —a~ 1 d. Thus, once a(s) is determined, 

jfe( s ) = _log^i and a(s) = a(s)- 1 (152) 
a 



in accordance with the general results of Section l2T6l 
4.3 Truncations of the hierarchy 

There are still infinitely many equations for the v m . We now study the accuracy 
of various truncations as a function of v 2 (ti) = A. Suppose we truncate by setting 



41 



v 3 = 0. Then v 2 (s) = A is a constant, and (|151|) integrates to 



as 



^/a§ + (x(0) 2 -x(s) 2 )4A . 
Similarly, the equations for A; and o, which in this truncation read 



fc 



4A !(*(*: 



2V 



a 2 , + (X(0) 2 - x(s) 2 )4A 



and <? 



4A X 



(a 2 , + (i 



x(s) 2 )4A ) 3 / 2: 



(153) 



(154) 



can easily be solved. Because x(0) = 1 and %(1) = 0, the final values ai = a(l), 
fci = k(l) and = o(l) are 



a i 



% + 4A , 



ki = log 



4Ar 



a 

^ = 72 



«0 



+ 4A, 



, (155) 



The thus obtained approximation to Po(h) is 



fcl+Sl 



|2_A 



1 + 



4A 



e a o +4A o 



l/,l 2 



( a o +4A o)' 



-H 4 



(156) 



By construction, this approximation is, for small Ao, exact to first order in Ao- But 
even at very large A 



e fe(i) 



v^+O(A 1/2 ) 



(157) 



— 1/2 

That is, the behaviour as a function of A comes out correctly as ~ A , only 
the prefactor is too small. Note that -Po(O) vanishes for large A , so the absolute 
error |Po(0) — e^ 1 ^ = 0(Aq 1//2 ) for large A . The difference of the true second 
moment to the one calculated in this truncation is much bigger: 

1 a 



wi - gi 



v/vfA^ ag + 4A 



(158) 



but the qualitative behaviour, namely the vanishing at large Ao, is still reproduced 
correctly. The fourth order moment has the right scaling as a function of Ao, but 
the coefficient is too large. 

The usual Polchinski hierarchy would not have given even a qualitatively cor- 
rect behaviour of these coefficients. Clearly, the adjustment of the propagator is 
crucial for large A because the end result W\ for the quadratic term is very differ- 
ent from the starting value 1 at large A . 
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However, the adjustment of the propagator alone is not sufficient for obtaining 
good asymptotics at large A because truncations at higher m improve the result 
only for small Ao, but not at large Ao- To see this, let us now truncate by setting 
v 4 = 0, so that (11481) becomes 



v 3 = Aa-'xvl (159) 

This corresponds to the usual truncation of the Polchinski equation, where the 
six-point function is replaced by the tree contribution (see Figure used to 
obtain the one-loop equation for the four-point function. In the general case, 
(|159|) is converted into an integral equation which is then substituted back into 
(|147|) . leading to a flow equation for v 2 that is nonlocal in the scale parameter s. 
In the toy model, ()147|) and ()159|) can be combined to 

9 

V 2^2 = --V3V3. (160) 

Because v 3 < at small s > by ()148j) . 



V3 = -\l^(X 3 -vl) (161) 



With this, (11471) becomes 



V2 = §a- 1 x^J-{\l-vl) 1 (162) 

and ()162|) and ()151|) now become a closed system. Now there is a serious problem: 
t> 2 decreases in the flow because % is negative, but this decrease never stops: even 
if t> 2 = 0, the A term in ()162|) keeps v 2 negative, and the flow then goes to 
negative v 2 and becomes unstable. A numerical analysis of higher truncations 
shows a partial stabilization, but there is still an instability already at moderately 
large A . The reason for this behaviour will be discussed in the next subsections. 



4.4 Discussion 

The above findings for the toy model can be understood easily in more general 
terms. The lowest truncation of the hierarchy is well-behaved in the toy model 
because the correction to a has a good sign, i.e. a(s) remains positive for all 
s because the A-dependent terms get added to it (in addition, it is reasonably 
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Figure 3: The six-point tree diagram gets integrated over scales s' < s and con- 
tracted with a propagator at scale s, to give the right hand side of the flow equation 
for the four-point function at scale s. 

accurate). The higher truncations become unstable at large couplings because 
the hierarchy is effectively nonlocal in the flow parameter: if one truncates by 
putting the eight-point function to zero, the right hand side of the flow equation 
for the six-point function t> 3 becomes a tree diagram of four-point functions, so 
that the equation can be integrated, and the equation for the four-point function 
v 2 becomes (assuming 1*3(0) = 0) 

JU 2 (s, *) = A D{s) j* ds' v 2 (s', *) F(s') v 2 (s', *) (163) 

The graphical representation of this procedure is shown in Figure |3] Because the 
entire flow history gets inserted on the right hand side, it does not help if v 2 de- 
creases in the flow because there is always the contribution from small s' where 
V2(s') is still close to Ao, hence large. This is the generalization of the right hand 
side of (| 162)1 . in which the Ao term kept driving the flow to negative v 2 even at 
v 2 = 0. Thus the nonlocality in the flow parameter makes it impossible to see 
in this truncation that the effective coupling constant decreases, but remains pos- 
itive, even at large initial coupling. In contrast, the 1PI equation is local in the 
flow parameter, and the just described "infrared asymptotic freedom" of the toy 
model at large A is seen easily. Thus in conclusion, to use the truncated equa- 
tions at strong coupling it is necessary that the scheme is stable in the sense that 
asymptotic freedom is not violated in these truncations even at large couplings. 

4.5 The Wick ordered flow in the toy model 

In the light of the above discussion, one may ask how an adaptive Wick ordered 
scheme would fare in this example, because the Wick ordered equation also has 
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the property that it is local in the flow parameter. As shown below, the Wick or- 
dered flow turns out to be well-behaved. In fact, already the initial Wick ordering 
provides a reasonable approximation to the true result at large Ao- 

As already discussed in previous sections, Wick ordering has to be done self- 
consistently if one wants to keep the interaction free of quadratic terms. The initial 
Wick ordering provides, when done self-consistently, the solution of the mean- 
field theory (Hartree-Fock equations) as a starting point for the flow. In the toy 
model, self-consistent Wick ordering means finding a such that 

aoM 2 + v{\(p\) = a\ip\ 2 + : v(\<p\) : a -i . (164) 

Because v(\ip\) = A |y?| 4 and : v : a -i= e~ a l9d v, the analogue of (j!08|) is 

|^| 4 = : ^l 4 : a _! +4a"V| 2 - 2a" 2 (165) 

Comparison of the quadratic terms in (j!64|) gives the self-consistency condition 

ao = a — 4a" 1 Ao- (166) 

The integral becomes 

P (h) = f^r e 2Aoa " 2 / dip A dip e -aM 2 -:AoM 4 : a -i+^+V (167 ) 
and the standard shift operations then give 

^(/f) = m e 2 ? e ^ fdfi a -i{(p) e- Xo ^+aW-^ (168) 



In contrast to mean-field equations in d > 1 dimensions, no integrations appear in 
the self-consistency relation (j!66|) . It is simply a quadratic equation with solution 



^iH 1 ^ 1 ' u69) 

For small A , (|169|) gives an O(A ) correction to a , and e 2A ° a 2 gives the lowest 
order correction to Po(0). For large A , 

a 



a = 2^X + ^- + ! J +O(X 1 ). (170) 



2 16^ 
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Originally, the strength of the coupling is given by the "dimensionless" ratio 
Ao/aQ. After the initial Wick ordering, this has decreased to 

£ = i-°«?) ("I) 

Thus the "strong coupling problem" in this integral is removed already by the ini- 
tial Wick ordering. Motivated by this, one can evaluate the remaining convolution 
integral in f!168|) to first order in A , to get the approximation 

p (/,)~M e 2Ao/a 2 e l^-^IM 4 (172) 

At large A , inserting ()170J) gives 

Po(0) ^7a=' w ^l x ^ /2 and < 173 ) 

Comparison with f!133|) and ()135|) shows that this approximation for \/Ao -fb(O) 

— 1/2 

differs from the exact value by 8 %, the one for the prefactor of A in w\ 
differs from the exact value by 12 % , and the approximation for the prefactor of 
A 1 in w 2 differs from the true value by less than 50 % . Note that these are the 
prefactors; the asymptotic behaviour as a function of A comes out correctly. Thus 
a simple initial Wick ordering, together with first order perturbation theory, gives 
a reasonable approximation for Po(h) even at arbitrarily large Ao- 

Let us now consider the simplest truncation of the Wick ordered flow with this 
initial condition. For the toy example, the choice r = d = aT x x 1S possible, so 
we use this to fix the Wick ordering covariance. It has the advantage that r is 
determined without any self-consistency argument, but on the other hand, some 
tadpole terms remain in the equation. Setting x(s) = 1 — s with < s < 1 
and expanding w(<p) = J2 m >2 w m\^\ 2m , the RG equations for the Wick ordered 
functions read 

a(s) = -(l-s)q(s) (174) 
itfa(s) = -20(1 - s)a(s)- 2 w 2 (sf (175) 
jfe(a) = -8(1 -s) 3 a(s)- A w 2 [\s) 2 (176) 
g(s) = a(s)~ 2 q(s) (177) 

The initial condition is given by the result of initial Wick ordering as (fixing a = 
1 in the definition of P and recalling ()169|) ) 

o(0) = a = ~ (l + v 7 ! + 16A ) 
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2A 

w 2 (0) = A , fc(0) = -lna+^f, #(0) = a -1 . 

and the function q(s) is given by 

( , 24 W2 (,) 2 a( g )- 3 (l - sf 
q{s) = l + 4a(s)-*(l-syw 2 (s) (178) 

The numerator is the value - in this toy model - of the two-loop sunset diagram 
which is the only one contributing to q(s) in the truncation w 3 = 0. The denomi- 
nator in the definition of q(s) comes from a tadpole term, which is there because 
d — f = —ar 2 ax- This is a special case of the term (1 — discussed in the 
previous section. Again, the sign is nice in the toy example — for u> 2 > 0, the 
denominator is always at least 1, so this inversion does not introduce singulari- 
ties. Only the equations (j!74|) and (|175|) are coupled; once their solution has been 
obtained, the other two simply follow by integration. 
Transform to the "dimensionless" functions 



«W = ^ and «(.) = ln^. (179) 

a(0) 



a(s) 2 



Then 



oj(s) = -lo(sY 20(1 - s) $[(1 - s) 2 uj(s)} (180) 

with = and w(0) = = \ - O(A x ). Given the solution of (ITSOl) . 
a(s) = a(0)e £ ^ s ^ is obtained from the solution of 

£(*) = -24^(,) 2 ? y (181) 

1 + 4(1 — syu(s) 

Equation ()180|) can be integrated explicitly but it is more instructive simply to 
note that cu(s) is decreasing and cu(s) > for all s because lu(0) > and the right 
hand side of ()180|) is negative and vanishes for u(s) — > (the last property failed 
for the Polchinski-type equation). Thus appears only with x E [0,o;(0)]. 
With w(0) = 1/4 this implies 0.7 < $(•) < 1, hence 

14(1 _ s )<_i^fl< 20(1 -s). (182) 

Thus the solution for u(s) is bounded above and below by that of the prototypical 
equation for "asymptotic freedom" . Solving this equation and setting s — 1 gives 
jj- < cj(1) < ^. Thus there is no instability in the Wick ordered scheme. 
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To obtain the coefficient of A 1 in w 2 , the correct starting values are cu(0) = 

1 /2 

1/4 and a(0) = 2A . With this, the truncated equations give the approximation 

W2 * a(OPeU) = 4^) A ° ' = °-° 33 A ° 1 (183) 



which differs from the exact value given in ()135jl by about 20 %. The solution for 
g gives W\ ~ g(l) = 0.578Aq Comparison with (I135j) shows a deviation of 
the coefficient of A 1 of about 2.5%. 

In summary, the bad behaviour of the Polchinski-type scheme at large cou- 
plings is absent from the Wick ordered adaptive flow. The standard truncation 
W3 = is stable and reproduces the coefficients of the expansion of the exact 
result in inverse powers of Ao well, i.e. with errors of order 1 - 10 % . Of course, 
the equation can also be used for all (not only large) values of Ao, where it has 
to be evaluated numerically. The stability properties are independent of A . The 
numerical values are in good agreement with the exact result. A similarly good 
strong coupling behaviour of the Wick ordered functions may also be expected in 
the more general examples discussed in [20]. 



4.6 Fermionic RG flows at strong coupling: general case 

In this section I briefly discuss the question of well-definedness of the RG flow 
in fermionic systems at large couplings. It has been stated in many places that the 
RG hierarchy is a diagrammatic method that cannot be taken as a starting point at 
strong coupling. The true situation is, however, not really so bad: the very strategy 
of the RG of integrating out degrees of freedom in small steps makes the RG flow 
well-defined also at large values of the initial coupling function. This is not a 
special feature of the differential equations but holds also for discrete RG flows, 
provided the stepsize is chosen small enough. 

To see this, consider a fermionic theory, as defined in the Introduction. By 
integrating the RGDE over short intervals, or using a discrete transformation, one 
can show that the generating function 

W(C, V)(0) = log J dfx c (<p) e xv ^ = log(/i * e" v )(0) (184) 

has the representation 

oo 

W(C, V) = nc * V + XP Wp( C , V, ... V) (185) 

p=2 
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where in W p , the argument V appears p times and 




(186) 

where T p is the set of all trees on p vertices, and P is a probability measure 
(i.e. positive and normalized, J dP(a) = 1), and M(o) is a positive definite 
p x p matrix with diagonal entries equal to 1. Both P and M can be writ- 
ten down explicitly, but only the above-mentioned properties matter (see ll29l . 
where the above formula is derived and explained in more detail). Eq. (I185|) 
can be thought of as a resummation of the perturbation expansion in terms of 
trees. Because ^\T P \ = — X ^P~ 2 < e p , the sum converges if the contribution of 
every tree T can be bounded by const p . For fermions, this can be shown us- 
ing determinant bounds (see, e.g. Il22ll29ll30l0 . The crucial point for the present 
discussion is that for a short-range interaction V sufficient conditions for conver- 
gence are that the determinant constant 5c of C (see El) and the decay constant 
otc = sup x J I C (X, Y) | dY are so small that | A| ac <C 1 . In a RG flow with step- 
size e, the covariance C is split into many small parts C £ and in every RG step, the 
C in the above formula is replaced by C £ . Because ac e ~ s and 8c e ~ s 1 ^ 2 , the 
above convergence condition can be satisfied for any A by making the RG stepsize 
e small enough. Thus, for small enough e the first steps in the RG integration are 
given by convergent perturbation theory in A, and therefore the RG flow can be 
started at large coupling. The really interesting analytical and physical question is 
whether a large, repulsive local initial interaction indeed decreases in the course 
of the flow (i.e. it gets screened), and whether the rate of this decrease is so fast 
that the flow can be controlled in the regime where the local coupling is not yet 
small, as is the case in the toy model. 

One may indeed expect that a strong on-site repulsion gets screened in the 
course of an RG flow, because even a hard-core condition gets softened, thus ef- 
fectively weaker when one considers larger blocks on the lattice. Namely, if there 
is a strong on-site repulsion, double occupancy of single sites is exponentially 
suppressed, but hopping into a block of m d lattice sites, m > 1, is still possible 
because of local fluctuations in the number density, even if the overall density 
is fixed. A standard example for a stochastic system with a hopping term and a 
hard-core repulsion is the asymmetric exclusion process. It was shown in [ 3 1 1 that 
averaging in a suitable way effectively allows to remove the hard-core constraint 
in this model and still get precise asymptotics of the decay of correlations. Be- 
cause the truncations of the correlation hierarchy in this example have the property 
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that they provide successive upper and lower bounds for the exact solution, one 
can obtain their asymptotic properties in some detail [ 3 1 1. For fermion systems, 
such monotonicity properties are in general not known, nevertheless averaging, 
such as provided by the RG, seems a very promising strategy to see screening 
effects. In fact, a low-order evaluation of the flow at very large scales in the two- 
dimensional Hubbard model shows the screening, as well as the generation of the 
usual antiferromagnetic interaction term 021 . 

The above argument is about the issue of the behaviour of an initial strong 
short-range interaction. The "flow to strong coupling" at low scales that has been 
observed in studies of initially weakly coupled systems is an independent phe- 
nomenon (as discussed in lfT3l . such flows may be kept under control for a while 
by phase space arguments, but eventually, the growth of the coupling function 
wins). The growth of certain parts of the coupling function in these flows is really 
due to the emergence of order parameter fluctuations that become more and more 
long-range, hence lead to fermionic interaction vertices that develop singularities 
in momentum space at a scale where symmetries get broken. In situations where, 
e.g. the singular part of the fermionic four-point function can be parametrized by 
an exchange of a boson that becomes massless at a certain scale, approximations 
that set the momentum of this boson to zero give a coupling constant of order 
m Boson which "runs to strong coupling" and diverges at a certain scale. Whether 
the boson can really become massless depends on the situation, in particular the 
dimensionality of the system. Moreover, in systems with long-range interactions, 
the Anderson-Higgs mechanism may interfere. However, unless this is the case, 
the Goldstone theorem implies such singularities whenever continuous symme- 
tries are broken. Because such singularities occur only at points, their effect is 
very different from that of a strong local interaction, which is equally large every- 
where in momentum space. 

5 Conclusions 

The RG equations developed here allow for a dynamical adjustment of propa- 
gators in a convenient and flexible way. In particular, they allow for countinu- 
ously adjusted scale decompositions, which are useful in Fermi surface studies, 
when one wants to avoid the use of counterterms, or fix the density (instead of the 
chemical potential) in RG flows. A numerical study of the adaptive flow equation 
discussed here, with fixed density, was done for the two-dimensional Hubbard 
model in the appendix of [8 1. A mathematical study for the case of regular Fermi 
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surfaces is in 0. Taking the Fermi surface deformation into account is important 
also for the case of singular Fermi surfaces, where the question of pinning of the 
Fermi surface at Van Hove singularities can be addressed using such flows. This 
situation is studied in Q using counterterms. The adaptive method is expected to 
lead to useful results also here. 

In flows with fermion and boson fields, a generalization of the strategy for the 
adaptive flow becomes obvious: quartic terms in the fermions contribute to the 
boson propagators, so the dynamical adjustment of the boson propagator allows 
to take out also certain terms from the fermionic four-point function. A related 
method has already been used in the 1PI scheme in ll33ll : a straightforward gener- 
alization of the method developed here will allow to do it in any scheme. 

Finally, the role of the propagator adjustment for the strong coupling be- 
haviour of flows was investigated for a very simple toy model. It was shown 
that the Wick ordered scheme with dynamic propagator adjustment correctly re- 
produces the asymptotic behaviour for arbitrarily large couplings, and already its 
simplest truncation gives the coefficients of the leading term in the asymptotic se- 
ries with an error of order 1 - 10 % (depending on the order of the moment). The 
analysis made clear that the adjustment of propagators is indeed crucial for the 
success of a scheme at strong coupling, but that in addition, a stability property 
of the four-point flow is necessary. Namely, the RG equations must be such that 
asymptotic freedom holds for the solution also if the flow is started at a large value 
of the coupling. The Polchinski scheme does not have the second property, even 
with an adapted propagator, but the 1PI and Wick ordered schemes do. 

Needless to say, the toy example is very special because there, the signs work 
out such that a large initial coupling Ao decreases to a very small final value, and 
the quadratic term ao increases to a large final value. Because only the combi- 
nation A/a 2 matters, this effectively leads to a very rapid transition to a small- 
coupling situation in the flow, provided the above-mentioned stability holds. This 
explains the success of the 1PI and Wick ordered scheme. In the latter, already 
the initial Wick ordering removes the strong coupling problem right away, if it 
is done self-consistently. It would be too much to hope for such benign signs in 
more general (and more interesting) models, but, as discussed in some detail in 
Section Ei~6l the flow equations are well-defined at strong coupling (contrary to 
what has often been stated in the literature) and there are reasons to be optimistic 
about screening of initially strong local repulsions in RG flows. 
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